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■ Abstract. This paper is devoted to the error estimate in the 
O^l , norm of the Adini element on the uniform mesh for the biharmonic 

equation. Surprisingly, a lower bound of the error estimate in the 
norm is established which proves that the convergence rate in 
the norm can not be higher than that in the energy norm. Let 
h be the meshsize of the triangulation, by the usual Poincare in- 

■ equality, this implies that the best error estimate in the norm 
I is at most of order 0{h^), which proves that convergence rates in 

both and norms are the same as that in the energy norm. 
This solves one long standing open problem whether convergence 
rates in and norms are higher than that in the energy norm, 
such a problem in fact dates back to the pioneering error analy- 
sis presented in [Lascaux and Lesaint, Some nonconforming finite 
[ elements for the plate bending problem, RAIRO Anal. Numer. 9 

\^ • (1975), pp. 9-53.] for this element. 
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1. Introduction 



For the numerical analysis of the finite element method for the fourth 
order elliptic problem, one fundamental problem unsolved is the error 
estimate in the norm [H [Sj [H]. In a recent paper [7], we analyzed 
^ i several mostly popular lower order elements: the Powell-Sabin — P2 

I macro element jl2j , the nonconforming Morley element fUX [T6l [TH [T7] , 

the — Q2 macro element [6], the nonconforming rectangle Morley 
element [15], and the nonconforming incomplete biquadratic element 
[T3l 118] . In particular, we proved that the possibly best error estimate 
in the norm for these elements is at most of second order and can 
not be two order higher than that in the energy norm. 

On the uniform mesh, the error estimate in the energy norm of the 
Adini element IH 131 [H [H] is of order 0{h?) provided that the exact 



Date: November 21, 2012. 

The research was supported by the NSFC Project 10971005. The research of 
the first author was supported by the NSFC Project 11271035, and in part by the 
NSFC Key Project 11031006. 

1 



2 



J. HU AND Z. C. SHI 



solution has the regularity [HI El [10] ; cf . [U [19] for the eigenvalue 
problem, where h denotes the meshsize of the triangulation under con- 
sideration. Consequently, the Adini element is in some sense a high 
order method for this case. Hence, by the usual Aubin-Nitsche trick, 
it seems that error estimates in and norms must be higher than 
0{h?). However, such an analysis can not be found in the literature so 
far. Therefore, it is still open whether error estimates in both and 
norms are higher than that in the energy norm for this element on 
the uniform mesh. 

The purpose of this paper is to analyze the error estimate in the 
norm for the Adini element PQ [21 [HI [Ilj • There are two main ingredients 
for the analysis. One is a refined property of the canonical interpolation 
operator, which will be proved by a new expansion method. The other 
is an identity for (— /, e)L2(Q), where / is the right-hand side function 
and e is the error. Such an identity separates the dominant term from 
the other higher order terms, which is the key to use the aforementioned 
refined property of the interpolation operator. Based on these factors, 
a lower bound of the error estimate in the norm is surprisingly 
established which proves that the best error estimate in the norm 
is at most of order 0{h?). Thus, by the usual Poincare inequality, this 
indicates that the best error estimate in the norm is also at most of 
order 0{h'^). Therefore, convergence rates in both and norms are 
the same as that in the energy norm and can not be higher any more. 
This solves one long standing open problem whether convergence rates 
in and norms are higher than that in the energy norm, such a 
problem in fact dates back to the pioneering error analysis presented 
in [Lascaux and Lesaint, Some nonconforming finite elements for the 
plate bending problem, RAIRO Anal. Numer. 9 (1975), pp. 9-53.] for 
this element. The result herein explains the numerical result in p]. 

The paper is organized as follows. In the following section, we shall 
present the biharmonic equation and the nonconforming Adini element. 
In Section 3, we define some interpolation operators and analyze their 
properties. In Section 4, we establish an identity for (— /, e)2,2(Q) and 
prove the main result that the error estimate in the norm has a 
lower bound which indicates that convergence rates in both and 
norms are at most of order 0{h'^). In Section 5, we end this paper by 
the conclusion and some comments. 
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2. The biharmonic equation and nonconforming Adini 

element method 

We consider the model fourth order elhptic problem: Given / G 
L\n) MweW := H^{n), such that 

(2.1) a{w, v) := {V^w, V'^v)^^^) = (/, v)L^n) for any v e W. 

where V^w denotes the Hessian matrix of the function w. 

To consider the discretization of (12. ip by the Adini element method, 
let 7/i be a regular uniform rectangular triangulation of the domain 

C in the two dimensions. Given ii" G 7^, let {xc, Vc) be the center 
of K, the horizontal length 2hj., the vertical length 2hy, which define 
the meshsize h := max(/ij., hy) and affine mapping: 

(2.2) ^ := ^^—j^, V ■= ^-j^ y) ^ ^- 

X y 

On element K, the shape function space of the Adini element reads 

(2.3) QAd{K) := PsiK) + span{xS, y^x} , 

here and throughout this paper, Pe{K) denotes the space of polynomi- 
als of degree < i over K. The nonconforming Adini element space Wh 
is then defined by 

Wh := {v G L'^{Q) : v\k G QAd{K) for each K ETh,v and Vf is continuous 
at the internal nodes, and vanishes at the boundary nodes on dfl^ . 

The finite element approximation of Problem f l2.ll) reads: Find Wh G 
Wh, such that 

(2.4) ah{wh,Vh) ■■= {Vlwh,Vlvh)mn) = if,Vh)L2(n) for any Vh^Wh, 

where the operator is the discrete counterpart of V^, which is de- 
fined element by element since the discrete space Wh is nonconforming. 

3. Interpolation operators and refined properties 

This section defines some interpolation operators and establishes the 
lower bound of ah{w — HhW, Htw) which is one main ingredient for the 
analysis of this paper. Here Uh is the canonical interpolation operator 
defined as in (13.31) below. 

Given K G Th, define the interpolation operator H^ : H^{K) — ?■ 

(3.1) {TIkv){P) = v{P) and V{Ukv){P) = Vv{P) 
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for any vertex P of K. The interpolation Ili^ has the following esti- 
mates [HEIEIE]: 

(3.2) \v - UKvlHiiK) < Ch^-'\v\H^iK)A = 1, 2, 3, 4 , 

provided that v G H*{K). Herein and throughout this paper, C de- 
notes a generic positive constant which is independent of the meshsize 
and may be different at different places. Then the global version Uh of 
the interpolation operator 11;^ is defined as 

(3.3) 11^1;^ = Ili^ for any K e Th- 
in the following, let V^f denote the i-th order tensor of all i-th 

order derivatives of v, for instance, i = 1 the gradient, and i = 2 
the Hessian matrix, and that are the piecewise counterparts of 
defined element by element. 

Given any element K, we follow O [7] to define Pkv G P^^K) by 

(3.4) / V^PKvdxdy= / V^vdxdy,i = 0,1,2,3,4, 
Jk Jk 

for any v G H'^{K). Note that the operator Pk is well-posed. It follows 
from the very definition of Pk in (13. 4p that 

(3.5) S/^Prv = Ho,i^VS 

with Ho,x the constant projection operator over K. Then the global 
version Hq of the interpolation operator Ho,k is defined as 

(3.6) Holi^ = Ho,x for any K e %■ 

Lemma 3.1. For any u G Pi{K) and v G QAd{K), there holds that 

3 Jj^ dx'^dy'^ dy"^ 

Proof. Let ^ and rj be defined as in (12. 2p . It follows from the definition 

of QAd{K) that 



(3.8) — = bo + bi^ + b2V + hir,. 



dy^ 
dxdy 



for some interpolation parameters Oj, bi, i = 0, ■ ■ ■ , 3, and Ci, i = 
0, ■ ■ ■ ,4. Since u G P^^K), the Taylor expansion and the definition of 
the operator yield 

n - Ilii-M = -(P - 1) + -T-TT-rl?? - 1) 
4! ; A\dy*^' ' 

hlhl d\ 



This leads to 



4! 9x4^ ^ ^ 2 dx^dy^^' 
(3.9) = 5^(12.2 _ 4) + ^J^U2 _ ^^ 

dxdy ^ ^ dx'^dy'^ 

A combination of (13.81) and (13.91) plus some elementary calculation give 

d^{u-IlKu)d^v , , hi f d^u d^v , , 
dxdy = — ^ I o.,9Q.,9 7r::^"^"^^ 



K 



dx"^ dx"^ 3 Jj^ dx^dy"^ dx 



(3.10) / 9'(" - n.n) |^ _ _| ^ 

ay^ oy^ 3 ox'^oy'^ oy^ 

which completes the proof. □ 

The above lemma can be used to prove the following crucial lower 
bound. 

Lemma 3.2. Suppose that w G H^{Q) fl H^{Q). Then, 

(3.11) ^h^ < (Vliw - Uhw), WlUhw)mn), 

for some positive constant (5 which is independent of the meshsize h 
provided that \'§;;^\h^{q.) 7^ '^'^d that the meshsize is small enough. 



(3.12) 
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Proof. Given K G 7^, let the interpolation operator Pk be defined as 
in (13.41) . which leads to the following decomposition 

= h + h. 

Let u = Prw and v = HkW in Lemma 13.11 The first term Ji on the 
right-hand side of (13.121) can be rewritten as 

3 Jk dx^dy^ dy^ 

This, error estimates of (13. 2p . and the commuting property of (13.51) 
yield 

3 Jk dx'^dy'^ dx'^ 

^1 f (^"^ ^"^^ d d 

3 Jf^ dx'^dy'^ dy"^ 

+ 0{h^)\\{I - Ti^)V^w\\L2i^n)\MH^n)- 

Since = (resp. = ) on the boundary parallel to the x- 
axis (resp. the y-axis), an elementwise integration by parts plus the 
uniformity of the mesh prove 



T — y n l|2 I x 

-'1 — ttW o 9o L2fo^ + — 



w 



(3.13) 3"ax29?/''^'(^^ 3"M?/2ll^'(f^) 
+ 0{h')\\{I -no)VM\L^^n)\\w\\HHn). 

The second term I2 on the right-hand side of (13.121) can be estimated 
by error estimates of (13. 2p and the commuting property of (13. Sp . which 
reads 

(3.14) I/2I < Chmi -Ilo)V''w\\L2^n)\\w\\HHn). 

Since the piecewise constant functions are dense in the space L'^{fl), 
II (/ - Uo)V^w\\L2(n) when h 0. 
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Hence, a combination of fl3.12p - p.14p proves the desired result. □ 

Remark 3.3. The expansion (13. 7p was analyzed in [ll[9l[T9]. Herein 
we give a new and much simpler proof. Moreover, compared with the 
regularity needed therein, the analysis herein only needs the regu- 
larity . 

Remark 3.4. The idea herein can he directly extended to the eigen- 
value problem investigated in [11 [19], which improves and simplifies the 
analysis therein and proves that the discrete eigenvalue produced by the 
Adini element is smaller than the exact one provided that the mesh- 
size is sufficiently small. In addition, such a generalization weakens 
the regularity condition from u G H^{Q) to u & if^(f2) where u is the 
eigenf unction. 

4. A LOWER BOUND OF THE ERROR ESTIMATE IN THE NORM 

This section proves the main result of this paper, namely, a lower 
bound of the error estimate in the norm. The main ingredients are 
the lower bound in Lemma [3^ and an identity for {—f,w — Wh)L'2{n) 
which will be established in the following lemma. 

Lemma 4.1. Let w and Wh be solutions of problems f l2.ip and f l2.4p . 

respectively. Then, 

{-f,w- Wh)mn) 

, =ah{w,^hW-Wh)-{f,^hW-Wh)L^isi) 
(4.1) 

+ ah{w - YlhW, w - Yihw) + ah{w - YihW, Wh - Uhw) 
+ 2(/, UhW - w)L2(n) + 2a/t(^i' - UhW, Uhw). 
Proof. We start with the following decomposition 

{-f,'w - Wh)L^n) 

= (-/, w - Wh) + ah{w, w -Wh) - ah{w, w - Wh) 

= (-/, w - Uhw)L'2(n) + (-/, ^hW - Wh)LHn) 

+ a{w, UhW - Wh) + ah{w, w - lihw) - ah{w, w -Wh). 

The last two terms on the right-hand side of (14.20 allow for a further 
decomposition: 

ah{w, w - Uhw) - ah{w, w - Wh) 

= ah{w - YihW, w - Tihw) + ah{TihW, w - Tihw) 

(4.3) - ah{w - YihW, w - Wh) - ah{IihW, w - Wh) 

= ah{w - UhW, Wh - Uhw) + ah{UhW, w - Uhw) 

- ah{UhW,w- Wh). 
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It follows from the discrete problem (12 ■4p and the continuous problem 
fl2.ip that the last term on the right-hand side of (14 ■3p can be divided 

as 

- ahOlhW,w- Wh) 

, = (/, ^hW - w)L2(n) - ah{w, UhW - w) 

(4.4) 

= (/, UhW - w)l2(s^) - ah{w - UhW, UhW - w) 
- ah{IlhW,IlhW - w). 

A summary of (I4.2p - (l4.4p completes the proof. □ 

Remark 4.2. The importance of the identity of (14.10 lies in that such 
a decomposition separates the dominant term 2ah{w — UhW ,Ilhw) from 
the other higher order terms, which is the key to employ Lemma \3.S[ 

For the analysis, we list two results from [8j and [9l [TO] . 

Lemma 4.3. Let w G Hq{Q) f] H'^{Q) be the solution of problem (12. ip . 
It holds that 

(4.5) \ah{w,Vh)-{f,Vh)\ < Ch'^\w\m{n)\\'^lvh\\L^n) for any Vh E Wh- 

Lemma 4.4. Let w and Wh be solutions of problems (12. ip and (12. 4p . 

respectively. Suppose that w G Hq{Q) fl H'^{Q). Then, 

(4.6) \\Vl{w~Wh)\\LHn)<Ch'\w\ 

Theorem 4.5. Letw and Wh be solutions of problems (12. ip and (12. 4p . 

respectively. Suppose that w G Hq{Q) fl H^{Q). Then, there exists a 
positive constant a which is independent of the meshsize such that 

(4.7) ah"^ < \\w - WhWmn), 

provided that \-Q^\m(n) 7^ and ||/||l2(q) 7^ and that the meshsize 
is small enough. 

Proof. We shall first bound the terms of the identity from Lemma [4.11 
The first two terms on the right-hand side of (14. ip can be bounded by 
Lemmas I4.3ti4.4l and estimates of (13. 2p . which leads to 

\ah{w, UhW - Wh) - if, UhW - Wh)mn) I 

< Ch'^\w\Hi(^Q)\\'^l(nhW - Wh)\\LHn) 

< Ch'{\\Vl{UhW - w)\\L^n) + \\Vl{w - Wh)\\mn)) 

< Ch^\w\l4^^y 



The estimates of the third and fifth terms on the right-hand side of 
(14. ip follow immediately from f l3.2p . which gives 

\ah{w- UhW, w - Uhw) + 2(/, HhW - w)L2(n) \ 

< Ch'^{\w\H4(^n) + \\f\\L2{n))\w\min)- 

From the Cauchy-Schwarz inequality, the triangle inequality, Lemma 
Ol and (O it follows that 



\ah{w - IlhW,Wh - Ilhw)\ < Ch^\w\ 



The last term on the right hand-side of (14. ip has already been analyzed 
in Lemma 13.21 which reads 



for some positive constant /3. A combination of these estimates states 

< {-f,w- Wh)mn), 

for some positive constant 6 which is independent of the meshsize pro- 
vided that the meshsize is small enough. This plus the definition of the 
norm of w — Wh proves 

{d,W-Wh) {-f,W-Wh) ,.. .. 2 

Ww-WhWmn) = sup —— > > 5/||/||L2{Q)/i . 

o^deL2(n) ll"l|L2(n) \\-j\\L2{n) 

Setting a = 5/||/||l2(q) completes the proof. □ 

Remark 4.6. For the rectangular domain Q under consideration, the 
condition \-§^\m(n) 7^ holds provided that ||/||L2(r2) 7^ 0. In fact, if 
\-§^\HHn) = 0,w IS of form 

w = coxy + h{x) + g{y), 

for some function h{x) with respect to x, and g{y) with respect to y. 
Then, the boundary condition concludes that both h{x) and g{y) are 
constant. Hence the boundary condition indicates w = 0, which con- 
tradicts with w ^ 0. 

Remark 4.7. By the Poincare inequality, it follows that 

ah^ < \\V{w - Wh)\\mn), 

for some positive constant a independent of the meshsize provided that 
|L2(n) 7^ and that the meshsize is small enough. 
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5. The conclusion and comments 

This paper presents the analysis of the error estimate of the Adini 
element in the LF' norm. It is proved that the best error estimate in the 
Li^ norm is at most of order 0{h^) which can not be improved in gen- 
eral. This result solves a long standing open problem in the literature: 
whether convergence rates in and norms are higher than that in 
the energy norm for the Adini element on the uniform mesh. This re- 
sult in fact indicates that the nonconforming Adini element space can 
not contain any conforming space with an appropriate approximate 
property. This will cause further difficulty for the a posteriori error 
analysis. In fact, the reliable and efficient a posteriori error estimate 
for this element is still missed in the literature, see [2] for more details. 
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